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Abstract 

An exactly solvable Quantum Field Theory (QFT) model of Lee-type is constructed to study 
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neutrino and the accompanying charged lepton can be calculated exactly as well as their creation 
probabilities. We can show that the coherent creation of neutrino flavor eigenstates follows from 
the common negligible contribution of neutrino masses to their creation probabilities. On the 
other hand, it is shown that it is not possible to associate a well defined "flavor" to coherent 
superpositions of charged leptons. 
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I. INTRODUCTION 



In recent times the neutrino oscillation phenomenon was established as an example where 
the coherent creation of a superposition of mass eigenstates, which is called a "flavor" 
eigenstate, tagged by the accompanying charged lepton, leads to an interference among 
those mass eigenstates inducing the flavor oscillation. We will choose to refer to those specific 
superpositions of mass eigenstates as neutrino flavor states because there is no well defined 
"flavor" charge operator to which we can associate flavor quantum numbers. Moreover, we 
will try to show that the concept of flavor for neutrinos is an approximate notion. 

In this article a Lee-tvpe model is deviced to investigate the neutrino flavor creation 
problem. The Lee model [l|, |2| is an example of a simple quantum field theory (QFT) which 
is exactly soluble but it still requires wave function and charge renormalization to make the 
model meaningful. Such features allied to strong conservation laws enable us to solve the 
model completely for each sector invariant by the conservation law. Another unique feature 
of this model, which will be explored here, is that it permits the study of unstable particles 
in an exact way 3j. The modifications performed here consist on replacing the fields N, 9, V 
by li, Uj, tt, where the unstable fermion V in the original model is replaced by the scalar 
boson tt and other fields are replaced accordingly. Thus this model will be used to mimic 
the decay ir ± — ► l^v a where a — e, /i. 

The decay process tt — > l^i/j, i,j = 1,2, has four kinematical decay channels available 
for two neutrino families. The channel tt — > l\v (li = e), with arbitrary neutrino content, 
is, however, very suppressed with respect to the dominant tt — ► l%v (I2 = /x) because of 
helicity suppression in weak decays. The model deviced here does not account for this 
helicity suppression, giving comparable branching ratios for the production of both charged 
leptons. However, comparable branching ratios for different charged leptons accompanied by 
neutrinos exist in nature l^: (a) K\ -> ^e^v e (40.53%) and K\ -> tt^^v^ (27.02%), (b) 
K+ -> 7r° e + v e (4.98%) and K+ -> 7rV + z^ (3.32%), and (c) W + -> l a v a . The last process 
is also effectively a two body decay kinematically analogous to the pion decay. 

In the aforementioned processes we could in principle observe the interference among 
neutrinos produced jointly with different charged leptons such as tt — > l\V\ and tt — > h^i- 
Why is this kind of interference not observed? We will try to address this question within this 
Lee-type model keeping in mind a fictitious pion decay without helicity suppression, with 
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comparable branching ratios to decay into different charged leptons. Previous discussions 
about the possibility of charged lepton flavor oscillations exist in the literature but 
we should remark that the question raised previously is more general in the sense that we 
could observe the interference of the channels by observing the neutrinos and not the charged 
leptons. On the other hand, there is also the interesting possibility that neutrino oscillations 
could be suppressed by small momentum or energy uncertainties (3, S]- 

A throughout study of the mixing phenomenon, either in the bosonic or fermionic case, 
was already performed in Quantum Field Theory (QFT) by Blasone and Vitiello (BV) 
Despite of that, the role played by the interaction which generates the mixing was not 
completely carried out. Although an analysis of flavor states emerging from pion decay is 



shown in Ref . 
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J. In an almost free QFT, BV analyzes the effects of mixing and seek for what 
should be the flavor states for neutrinos. What emerges from that study is the possibility to 
define two types of unitarily inequivalent states constructed from two distinct vacua: a usual 
vacuum and a mixed vacuum. They use the mixed vacuum and obtain a different neutrino 
oscillation phenomenon. The choice of the vacuum, however, does not seem to be unique. 
The study of a Lee type model presents two advantages with respect to the BV approach: 
(i) the interaction responsible by flavor mixing can be taken into account exactly and (ii) the 
free states are unitarily equivalent to exact states of the total Hamiltonian. The point (ii) 
means we do not have to choose between two unitarily inequivalent set of states but at the 
same time we can not study in this type of models the consequences of such phenomenon, 
which might be present in more realistic QFTs. 

Another distinctive approach to treat neutrino oscillations in a rigorous manner was ex- 
tensively reviewed in Ref. [ill, where the so called external wave packet (EWP) approach was 
analyzed. Although, the first QFT treatment of neutrino oscillations was given in Ref. [12. 
In the EWP approach, the nonobservability of neutrinos through direct detection was trans- 
lated into the formalism by treating them as internal Feynman propagators connecting the 
creation and detection processes. Some unclear aspects of such approach was clarified in 



Ref. 



13l . One of them concerned the possibility that some nonphysical contributions, such as 



below threshold neutrino detection, could be present in the formalism and they have to be 
removed by hand j3| . Such approach could account for the localization aspects of the cre- 



ation and detection processes 



14 



15 



16] but the meaning of what would be the intermediate 



neutrino flavor states could not be investigated. Another attempt to define neutrino flavor 
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states in QFT can be found in Ref.ll7l. One can also find attempts to describe neutrino 
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oscillations in Relativistic Quantum Mechanics 

The two goals of this article consist in providing an exactly solvable model to study (a) 
how neutrinos are created through interactions and (b) how the localization of neutrinos 
follows from the parent pion properties. As is usually adopted in rough estimates the wave 
packet size of the daughter particles is considered to be 1/r [5|, [7|, where T is the decay width 
of the parent particle. This model provides exactly such relation. We can also calculate the 
two particle state composed by the charged lepton and the neutrino that is created from 
pion decay. This calculation is an attempt to properly define a flavor state for neutrinos. We 
will see that the coherent creation of neutrino flavor states is a consequence of the common 
negligible kinematical contribution of the different neutrino mass eigenstates to their creation 
probabilities (amplitudes). We will also study the meaning of defining a different "flavor" 
state to charged leptons as the coherent superposition that accompanies a neutrino mass 
eigenstate such as u±, for instance. We will conclude that the concept of "flavor" associated 
to such superposition states can not be properly defined. 

The outline of the article is as follows: in Sec. [Ill we introduce the model and calculate 
the relevant eigenstates of the total Hamiltonian. In Sec. lIIII we study how the concept 
of neutrino flavor emerges from the calculations and derive the localization properties of 
the decaying states. The discussions are made in Sec. [IV] and some detailed or additional 
calculations are shown in the appendices. 

II. THE MODEL 

The free Hamiltonian of the model is defined by 

H = (i J d%At(k)A(k)+M i J d^6t(p)6 j (p) 

+ J d^^(p)ot(p) ai (p) , (1) 

where 

E j (p) = y Jj? + m}, (2) 

and A, hi and cij are respectively the annihilation operators of the fields it, U and Uj modelling 
(unrealistically) the decay it — > The summation over repeated indices are implicit. 

We will restrict ourserves to two families which means U = h,h {h = e, I2 = fi) and 



Uj = V\,V2 correspond to the two neutrino mass eigenstates. Extension to three neutrino 
families is straightforward. The fixed bare energies of the fields n and Zj corresponding to 
their masses fi, Mi mean this model treats the static and recoilless limit of tt and U. The 
momentum though is conserved. The neutrinos, on the other hand, obey the relativistic 
energy dispersion relation (J2J). The creation and annihilation operators satisfy 

[A(k),At(k')] = 5 3 (k-k') (3) 
{h{p)Mp')} = M 3 (p- P ') (4) 
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The spin indices are suppressed. 

We define the Fourier transform of the fields to be 



{a,(p),a](p')} = M 3 (P"P) (5) 



e 

k{x) = I d 3 pbi(p,r)u r — -y (6) 



_.J2 

T(X) = / AA(k)-^ (7) 



(2k) 



—=[a i (p,rX(p)-—j + al{p,r)rjcCu r i (p)— — T ] , (8) 

V^A (27TJ2 (27r)2 

where «[(p), r — 1,2, are Dirac spinors normalized as u['(p)«|(p) = 2i?j(p)5 rs while Mq, 
r = 1,2, are the corresponding spinors for fermions at rest properly normalized to UqUq = 
8 rs . Expression (jHJ) corresponds to the usual Majorana type fermion expansion where rjc is 
a phase factor appearing in the charge conjugation transformation. 
The interaction Hamiltonian is chosen to be 

Hi = g J d 3 xd 3 y [4(x)^z/j- ) (x)7r(y)/(x - y) + h.c] (9) 

= goUij J d 3 pd 3 p'b\(p, r)$V)aJ(p'> *)Ml> + p') + ^ ( 10 ) 

where {Uij} corresponds to the mixing matrix, z/j _> refers to the creation part of expansion 
© and 

CV) = * a i^ ^o^V) • (11) 
(2tt)2 yj2E j (p>) 

Using the Dirac representation for gamma matrices, C = 27270, Uq = e r [19J] and Uj(p) = 

E„^ +mj ^ obta . n 

n,,(P) = n^Ci? T (p) = -'"^'f;'" . (12) 
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The function / represents a "form factor" necessary to regularize the expressions in the Lee 
model j^j]. For consistency reasons to be discussed later it also has to be smooth, analytic, 
approximately flat for energies much lower than the cutoff scale, but falling off not so rapidly 
above the cutoff scale. In addition, the Lee model may violate unitarity if the coupling 
constant is above a critical value • We assume the coupling constant is below the critical 

ess, the case of an unstable parent particle can 



3] 



value thus avoiding such situation. Neverthe . 
be properly described in the unitary regime 

There are two charges Qi, Q 2 that are conserved and within those sectors of fixed charges 
we can find the explicit eigenstates of the whole Hamiltonian 



H = H + H 1 . (13) 

The two charges are 

Qi = J d 3 kA^(k)A(k) + J2 f d 3 pbl(p)h(p) } (14) 

Q 2 = J2 I <5»[6j(p)6i(p) - 4(pVi(p)] • (15) 

In the sector containing one n [(Qi, Q2) = (1, 0)] and one pair of U, Vj [(Qi, Q2) = (1, 0)] 
we can calculate the exact eigenstates of H in terms of the following eigenstates of Ho 

|tt) ee J d 3 A:^(k)|7r(k)) , (16) 

\k(r),Uj(p',s)) ee J d 3 pi) 7r (p + p')\l i (p,r),v j ( y p,s)) , (17) 



where 



|7r(k)) ee A+(k)i0>, (18) 
|/ i (p,r),^(p', S )) = 6l(p,r)at(p',s)|0), (19) 

and |0) is the vacuum state of Hq and of H as well [20]. The function ip w is an arbitrary 
function controlling the pion momentum distribution. The subscript indicates the states 
correspond to eigenstates of H , i.e., bare states. Then we find the eigenvalue equations 

#ik)o = goUij J dVtf/ip'Mr), uj{p', s)) , (20) 
F 1 |/ J (r),^(p', S )) = ^ ^r(P')k)o- (21) 
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The eigenstate of H with energy (mass) E w corresponding to a dressed state of \tt) can 
be found if we discover the function x in 

k) = k>o + X) / dVx^(p')IW,^(p',a)>o- (22) 

From the eigenvalue equation 

H\n) = E w \7r) , (23) 



we obtain 

x -(p') = g^gg^ (24 ) 

and ho(E n ) = where 

^(K) = E % - // + O (K) ■ (25) 

The function O is defined as 

rs 



We assume there is only one root for Eq. (125]) which defines the pion mass M n , i.e., the 
energy of the dressed state \tt). Therefore we could have written the bare pion mass fi as 



where 

is the mass counterterm. 
We can simplify 

where 



fx = M n + Sfx, (27) 



5 fx = MM,) (28) 



^(E)= l0 Y,\U ij \ 2 4> j (E-M i ), (29) 



7o^ (30) 



and 



(f)j{x)= I dE x lE 2 -m)(E-m 3 )Q^-. (31) 



The intermediate steps to reach Eq. (1291) can be found in appendix El We should emphasize, 
however, that the principal value is understood in Eq. (!3~T|) if x > rrij j^]. 
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To have an unstable pion we 



should have |3 



h (M 1 +mi) < 



(32) 



for Mi = min(Mj) and mi = min(m,,). This condition implies that 



M w > Mi+mi. 



(33) 



These conclusions follow from <p (—oo) = — e, d(j) (x)/dx > for x < mi + M x and 
4>o(x) > for x — > mi + Mi — e, which assures h^{x) is a monotonically increasing function 
for x < Mi + mi. Equation (1321) guarantees that the root of ho(E) = 0, the pion mass, is 
above the threshold Mi + mi, leading thus to Eq. ( 1331) . It can be also seen from Eq. ( 13T1) 
that for x > Mi + mi, 0o( a; ) ma Y continue increasing but it starts to decrease when x is 
greater than the cutoff scale of /. In addition we will assume all the decaying channels are 
open, i.e., M n > m 2 + M 2 . Intermediate cases can be handled with appropriate care. 

If condition ( J33l) is satisfied, Eq. (!22l does not define a meaningful expansion although 
the norm is finite 



with h' Q (x) = dh (x)/dx. For stable pion, the norm (|3"^1) is used to (re)normalize the pion 
state since the expression in Eq. diverges when there is no cutoff, i.e., f(E) = 1. It 
means there is no stable dressed pion state and the only stable eigenstates of H in this sector 
(one 7r or one pair of li,Vj) are the scattering states containing one U and one Vj. These 
states complete the Hilbert space in this sector. This is proved in appendix IB"1 
The scattering states, eigenstates of H, can be calculated from the expansion 



(34) 




i 



+ Z*03(p)\ir) o , 



(35) 



if we find the functions a and Z| (3. The eigenvalue equation 



H\k{r), Uj(p, s)) = (Mi + Ejip)) \k(r), ^(p, s)) , 



(36) 



leads to 



(M f + E,(p) - - ^,(p'))ay,5'(p, p') = gcZ}^(p)4 , '(p')U i , j , , (37) 
(Mi + Ej(p) W^(P) - «HPW + * / (38) 
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Inverting Eq. (137)1 using the incoming (+) or outgoing (-) wave boundary conditions we 
obtain 

Equation (1381) yields 



Zj^(p) = »£» . . (40) 

Notice that the domain of the function ho was extended to the complex numbers, hence 
satisfying the property 

h {E±ie)=h {E)± l -V (E), (41) 

where 

T (E) = 7o ^)|^| 2 f i (E - Mj) , (42) 



fj(a;) ee 2n6(x — rrij)(x — irij)Jx 2 — rrvjf (x) . (43) 



Therefore 



a -v (±)( , = ggM! ^ (44 ) 

Despite the divergence in Eq. (13~4"|) . we can still define the approximate pion state 

r ri rs (v') 
M=C[ W „- Sot /« j dy Mi + ^: M ^_ iX Hr)Mv',sM, (45) 

with M n being the root of Eq. (|25p . Such state can be properly normalized if A 7^ 0. That 
this state is not an exact eigenstate of H is proved in appendix [Dl If we expand the state 
( l45l) in terms of the scattering states ( 1331) . we can see how this state can be seen as an 
approximate ressonant state corresponding to the pion that decays. 
We calculate, for outgoing states, with superscript (— ) suppressed, 

(k(r),u 3 (p,s)\n) = Z*l%"(p) (46) 
(i»(r),Vj(p,s)|/i/(r),i/j/(p',s))o = <5 3 (p - p')drr'S 8a >5ii'8jj> , (47) 

which yield 

goUijVj'iP) 



rs 
ij 



h {Mi + Ej + ie)(Mi + E j — M n — iX) 



X [ - l\ + 0o (Mr) " 00 (M w + iA)] . 
For small A, i.e., |A| <^ r o (M 7r )/20Q(M 7r ), we can approximate 

r(M 7r ) 1 l 



(48) 



|tt a ) « -iA 



1 + 



2IAI 



rs 



E / d *P\ l i(r),Vj(p,s)) 



h(Mi + Ej + ie)(Mi + E 3 — M n — iX) 

We have denned the renormalized coupling constant 

i 

9 = Z%9o, or 7 = Z 2 7o, 



(49) 



(50) 



which defines the "renormalized" functions I, obtained by replacing 7 by 7 in T , 0o, thus 
multiplicatively, while 

h(E) = Z 2 h (E) . (51) 

The renormalization constant is 

1 



Z 2 - 1 = l + 0[ ) (M 7r ) 



(52) 



1 - / (M 7r ) ' 

Since the scattering states \h(r), s)) are eigenstates of the total Hamiltonian H with 
energy E = Mi + Ej(p), we can easily compute 



1 + 



r(M 7r ; 
2|A| 



mi+Mi 



2vr 

-iEt 



[(£ - + ^\E)f + JIM] [(£ - M^f + A 2 ] 



(53) 



where |7T A (t)> = e~ iHt \Tix), (2) (£) = <p(E) - 0(M W ) - (j)'{M^)(E - M w ) and Eq. were 
used. 

If the functions T(E), <p(E) vary more slowly than E near = M n we can approximate 
T(E) ~ r(M vr ) and if M„- ^> Mj + we can replace the lower limit of the integral in 
Eq. (E3D by -00: 



(irx\irx(t)) « e" 

|A|>r w 



1 + 



r, 



-1 



where 



2|A| 

-i(M- IT -iT 7T /2)t 



T n = T(M V ) 



-rvt/2 _ r 7r -\x\t 
2IAI 



(54) 
(55) 

(56) 
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Equation (1531) expresses the usual exponential decay law with decay rate IV This result is an 
indication that the approximate state \tt\) in Eq. (j4"5"l) represents appropriately a ressonant 
pion state that decays into \h{r), z^(p, s)) states independently of the arbitrary parameter A. 
If we had calculated Eq. ( 1541 for a state analogous to Eq. ( )45i) but with approximate energy 
very different from M„ we would obtain a decaying amplitude in Eq. ( I55I) with decay rate 
primarily given by |A| ^> T n (see proper discussion in Ref.[3|). It is important to remark that 
|A| should be also sufficiently small to allow Eq. 04 II) when e is replaced by A. Combining 
|A| > T w and |A| < T w /2(f)'(M n ), it is necessary that (f)'{M^) < 1. 
The emission spectrum of k , Uj is given by 

|(/ J (r),^(p, S )|7r A (t))| 2 =|C17(p,A)| 2 , (57) 

where 



gUijVj a (p) 



C-(p, A) = -iX 1 + — ^ , - , - - ^ (58) 



h(Mi + Ej + ie)(Mi + E* — M n — iX) 



is the expansion coefficient in Eq. ([49]) . We can see 

/poo 
d*p\C™(p,X)\ 2 = J2 / dEW lJ (E-M l ,X) = l. (59) 

ij 



The function Wij(E, A) is the emission probability for n — > Uvj, per unit energy, with 
neutrino energy E, summed over spins: 



l2 Wx(E) 



W l3 (E - M h A) = r yTj(E - Mi)\U^ . (60) 

The total emission probability is 

W\(E) = -1 ^ A2 [ 1 + r -/ 2 ' A 'J (61) 

XV ' 2tt h 2 {E) + T 2 {E)/ 4 {E-M n ) 2 + X 2 1 ' 

which approximates the usual Breit-Wigner distribution in the A ^> T n limit, within the 
approximations h(E) ~ E — M n and T(E) ~ r„-, 

^).lim^( g ) = ± (g _ M ^ 2 + r;/4 . (62) 

From probability conservation we have 

Wij(E - Mi, A) = W X (E) . (63) 

The energy E in W^i?) refers to the total energy of the pair I + v. 
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We can see the pion state is a ressonant state by calculating the scattering cross section 
of l; L Vj — > lu, irrespective of the final product flavors and spins and averaged over initial 
spins: 

-* lv ' E) = (E-S-^ W ~ M '» i w{E) • (64) 

where E = Ej + Mi corresponds to the total energy of the incident particles. The kinematics 
is totally determined because the energy of k and tc are fixed. The difference between the 
initial and final momenta (velocities) of neutrinos has to be taken into account to obtain 
the appropriate cross section obeying 



p2=( jB _M l ) 2 -mf 



The presence of W(E) guarantees the resonance for E = M n . 

III. FLAVOR (EIGEN) STATES 

In the previous section we have seen, from the QFT point of view, it makes more sense 
to calculate transition probabilities (amplitudes) with respect to the state \li(r),Uj(p,s)) 
that has definite energy Ej(p) + Mj and thus are mass eigenstates. What wee see in the 
actual pion decay, however, is not tt — > 1(1/ j, i,j = 1, 2, but n — > l a v a , with a definite flavor 
a — e, fx. For the charged leptons the terminology "flavor" is clear and coincides with the 
mass eigenstates U = l a . For neutrinos what characterizes a flavor state depends on the 
accompanying charged lepton. 

Let us define flavor states as 

\l a (r),up(p,s)) = SatUpjlkW^jfas)) , a,p = e,fx. (66) 
Let us then calculate 

\(l a ,vp(p)Mt))\ 2 = ^|^^^(P)^ (Ma+£j(p))t | 2 (67) 

rs j 

=»'EiE^'(^ MfiJ T£ +fe) r. <« 8 > 

already in the \ ^> T w limit. If we take the rrij — > limit (Ej — > E v = |p|) in the terms 
except in the exponent we arrive at 

\(L,Mp)\n x m 2 = V a ,(t)^^, (69) 
12 



where 

V af 3(t) = \j2Ua j e- iE * t U; j \\ (70) 

3 

and W a (E) is W a j(E) in the rrij — > limit and with |t4/| 2 factored out. We see the 
probability of finding the flavor state (l6"oT) from pion decay decouples into the usual oscillation 
probability times the emission probability it — > l a v irrespective of the neutrino flavor. 

The concept of neutrino flavor states emerges from Eq. (1691) taking t = (or t <C oscillation 
period), for which we obtain 

V afS (t^0)=5 af3 . (71) 

Therefore only a neutrino of flavor a is produced jointly with the charged lepton of flavor 
a. If we try to calculate the creation probability for the equally possible superposition 
state, corresponding to a neutrino mass eigenstate accompanied by a superposition state of 
charged leptons, 

\J2 U h( l ^3(p)Mt))\\ (72) 

i 

we see it does not decouple and since the masses Mj have different kinematical contributions, 
there is no way to define a superposition state of charged leptons associated to a neutrino 
mass eigenstate vi for instance. Stated differently, the transition probability ( IT2l) is not 
equal to 8pj for t = 0, not even approximately. 

The definition of the neutrino flavor state in Eq. (1661) was justified only a posteriori as 
the state that accompanies a definite charged lepton. In this model, however, we could, in 
principle, calculate the more interesting quantity 

(^(rV^x^Mi)) = %^t)e- m >\ (73) 

where 

\li(r),Vj(x,s))= I d 3 p^— T |Zi(r),^-(p,s)). (74) 



(2tt 

This "wave function" squared (|\l/^(x, t)| 2 ) would give us the probability density of finding 
the neutrino Vj in the position x at time t jointly with a charged lepton Zj. The exact 
information of \P£?(x, £) allow us to expand 

M*)) = E/ d\%^t)e-^%{r)^ {^s)) . (75) 

ij,rs 
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Equation (1751) characterizes completely how the neutrinos and charged leptons arise from 
pion decay. Although the spatial information of the charged lepton is not determined or 
calculable in this model. 

The exact Fourier transform necessary to compute Eq. (173]) is very difficult to be per- 
formed analytically but a rough approximate calculation can be performed. We find, in the 
limit |A| ^> rV, 

= Jd'pe-^C^ip)^, (76) 
' ~' L '' 7 i^f)(4-) ( °"^ 2)rS \6(r-)e^ r e-^ T - + 6(r + )e-^ r e-^ T +] ,(77) 



where 



r ± = t± r/vj , (78) 

r = |x| and the bar in Vj,pj and Ej denote respectively the velocity, momentum and energy 
of Vj when Ej = M n — Mi, which corresponds to the kinematics of it — > liVj with tt, li at 
rest. The detailed calculation of Eq. (1771) is carried out in appendix O 
We can verify the normalization of Eq. (1771) : 

[ d 3 xm°(x,t)\ 2 = ^|^ff^^.)[l + 4^e- r ^sin(^)]. (79) 
Equation ( !79l) guarantees that 

£ [ d 3 xW;(x,t)\ 2 = l, (80) 

ij,rs 

for t ^> l/TV or ^ <C 1. Therefore, the total probability is conserved by the approximation. 
One can not, however, recover the flavor oscillation phenomenon from the approximate wave 
function ( 1771) because the phase is common to all states Uvj with fixed i. Numerical studies 
can be performed to test the accuracy of the approximation, although appropriate care 
should be taken with the enormous difference between the scales of rrij and M n , Mi. 

To visualize the localization aspects of the neutrinos created, we can calculate the radial 
probability density for the creation of neutrinos Uj at radius r and time t jointly with l iy 

Pij (r,t) = W dQ x r 2 |^(x,t)|, (81) 

rs 

= 7 |^| 2 f,(M 7r -M,)-^(r_)e- r ^- 



V 3 
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+ e(r + )e~ r ^ T+ + 26(t + )6(t_) cos(p i r)e- IV *] . (82) 

For t ^> 1/T n , we see only the first term containing r_ contributes. Thus Pij(r,t) has a 
triangular shape with an abrupt peak in r = Vjt and a exponential tail for r < u,t, being 
negligible for r > Vjt. The neutrino is then roughly localized in the region Vj(t — 1/T n ) < 
r < Vjt. The size of the wave packet is roughly Vj/T n as is usually assumed in rough 
estimates js, 3|- 

The most faithful state that we can construct to describe the neutrino flavor states created 
from pion decay jointly with a charged lepton Zj with momentum q and spin r is 

AMit 



*o<«i(q,r)|7r(t)> = E/ ^(q + pK^ (p) *Mp, s)> 

gUijr)j S (p) 



3,8 

x 



e -i(Af T -iWi-£?j)t e -r w t/2 



83) 

^(Mi + E^p) +ie) r " J v 

where ^/v is the function appearing in the definitions (116ft and (ITTjl . corresponding to the 
pion momentum wave function. In the limit t ^> F % , we can approximate rrij — > in all 
factors except in the exponent since rrij <C |p|,Mj, M n . Equation fl83|) becomes 

e m \U^r)\ix{t)) nj^j d^M + p)e-^ t U lJ \u 1 (p, a))o fe(M ^ ~ ) • (84) 
We could then write, for £ = 0, 

o<Ja(q,r)|7r(f))«X) | dW»(q + p)|^(p, s)) j^f^ \ - e) , (85) 
where h = l a , Mi = M a and 

k«(p, s)) = Y U <tf\ u j('P> s ))o • (86) 

3 

If we impose Vv(p) = <5 3 (p) we obtain exactly q = — p for the momenta of the charged 
leptons (q) and neutrinos (p). In general, the presence of the function if> n ensures momentum 
conservation and forces the amplitude fl83l) to be appreciable only around q = — p where p 
has magnitude satisfying Ej(jp) ~ M n — Mj. Equation (1831) also clarifies the contribution 
of the intrinsic momentum uncertainty Ap (intrinsic of of the parent particle that is 
inherited by the daughter particles, apart from the contribution of the decay width rV the 
smallest between r„- and Ap dominates. Other transition amplitudes of \ir(t)) with respect 
to the free states \tt)q and \k(r), Vj(p, s))o can be found in appendix [El One can identify 
some similarity with the Wigner-Weisskopf approximation used in systems with couplings 
between discrete and continuum energy levels. 
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IV. DISCUSSIONS 



One knows neutrino flavor oscillation ceases when quantum coherence is lost due to the 
lack of spatial overlap among the mass eigenstates that compounds the flavor state neutrinos 



21 



22] . The characteristic time (distance) scale for such phenomenon to occur is usually 
very large for neutrinos and is given by Sx/Avij, where 6x is the characteristic spatial size of 
the neutrino wave packets and Avij is the velocity difference of the Vi and Uj that compound 
the flavor state. In this model, in the pion restframe, 

5xj « Vj/T n . (87) 

Neutrino flavor oscillation occurs because the neutrinos produced in channels tt — > l 2 v\ and 
7r — ► l 2 v 2 interfere coherently and remain interfering in the same region of space as long as 
t <C 5xj/Avij. But what prevents the channels tt — > l\V\ and tt — > l 2 u% from interfering? 
For real pion decays such interference is not observable even in principle since these two 
channels have very different probabilities to occur because of helicity suppression in weak 
decays. In this model, however, the branching ratios are comparable and may mimic other 
real decays. 

One should notice that the neutrinos created in channels tt — > l\V\ and tt — > l 2 V\, for 
example, may have significant spatial overlap. We see that AMv(rrij) defined as 

A M v( mi ) = v^Mt) - ^(M 2 ) , (88) 

may be of order of 

A m ^(M,) = Vi(Mi) - v 2 (Mi) , (89) 

unless the neutrino masses are nearly degenerate Am| <C m^,m^. We use the notation 
Vj(Mi) instead of Vj throughout this discussion to avoid ambiguities. It means that the non 
interference of these channels is probably due to the charged leptons, but its detailed account 
deserves further study. Simple calculations show [g] that the big mass difference between 
charged leptons make them lose spatial coherence over a distance of atomic length or less for 
most decaying processes. Is this lost of spatial coherence responsible for the incoherence of 
the channels with definite flavor li — e or l 2 — /x? Another possibility is that the detection 
of charged leptons is what makes the coherent superposition of the four channels it — > l^Uj, 
i,j = 1,2, reduce into two incoherent n — ► l\v and tt — > l 2 v quantum processes. These two 
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explanations are, however, very distinct since the former is detection independent while the 
latter is detection dependent. 

Concerning the absence of the interference between the channels 7T — > l\V and tc — > fav, 
there is one ingredient that could play an important role, i.e., the spatial entanglement. 
This question, however, can not be discussed within this model due to the lack of spatial 
entanglement between Zj, in the static limit, and Uj. 

Another known situation when neutrino oscillations can be suppressed is the case when it 
is possible to probe the neutrino kinematical quantities with enough precision to single out 
one neutrino mass eigenstate that is being created in the process [71] • In the language of this 
article it means all the four pion decay channels proceeds incoherently. Such suppression, of 
course, can be induce d by localization properties of the source or detector that produces and 



detects the neutrinos 



141 ] . Conditions necessary for neutrino oscillation can be analyzed, for 



instance, within a realistic OFT description 16}]. In that case, the condition for neutrino 

n 

oscillation requires [Eq. (4.7) of Ref. [16[] 

if the intrinsic momentum spread of the decaying particle, apart from the decay width F, 
is small compared to its mass. Condition (1901) is exactly the condition necessary to prevent 
one from knowing which neutrino mass eigenstate is being created [7j], when the momentum 
uncertainty is of the order of T, which is also the case of this work. Notice, however, that 
in this work we do not treat the propagation of neutrinos but rather focus the attention on 
the creation process. Surprisingly, condition flUUl) . with <C instead of <, is exactly what is 
needed to guarantee, initially, approximate neutrino flavor definition. Such analysis involves 
estimating the first flavor violating term of Eq. f l68l) that follows after Eq. ( |69i) . 

Therefore, the study of this Lee-type model sheds light into the question of what is a su- 
perposition state of definite flavor. We could confirm, within this model, the result that once 
the channels with charged leptons with definite masses are singled out, coherent creation of 
neutrino flavor states is allowed because the kinematical contributions from each neutrino 
mass eigenstate to their creation probabilities are negligible. Therefore the creation proba- 
bility amplitudes of neutrinos are mostly independent of neutrino masses (mass differences 
to be precise), enabling the decoupling of Eq. ( |69i) into the creation probability times the 
flavor conversion probability. The conversion probability, however, depends crucially on the 



17 



tiny neutrino mass differences. The neutrino ffavor states u a can be created jointly with the 
charged lepton Zj = l a and the wrong ffavor is initially absent. On the other hand, if we try 
to calculate, as in Eq. (1721) . the creation probability of the coherent superposition of charged 
leptons jointly with the neutrino v\, for instance, there is no meaning to attribute the name 
"ffavor" for that superposition of charged leptons because such ffavor is not created in an 
unambiguous manner but both the correct and the wrong flavors would be created initially. 

To be extremely precise, the notion of neutrino flavor states rely on the approximation of 
negligible neutrino masses, as could be seen in Eq. (1691) . For that reason, there is a negligible 
probability to create the wrong flavor v ei for example, jointly with the charged lepton I2 = fi. 
Such probability is calculable but negligible because of the proportionality to the neutrino 
masses. Such flavor indefiniteness due to the neutrino mass differences can be also seen in 
other attempts to exactly define what is a neutrino flavor state such as in Refs.y or 17. 



The exact definition is always blurred by the differences in the kinematical or dynamical 
factors that depend, strictly speaking, on the neutrino masses, that in turn, can not be 
factored out as a common factor. For example, the spinorial character of neutrinos prevent 
even a description as simple as a mixed free second quantized theory to provide exact initial 



neutrino flavor definition 



131 ] . Another simple treatment that leads to initial flavor violation 



is the adoption of general kinematical conditions when scalar wave packets are used 23]. 

To summarize, we could confirm two aspects, within a Lee-type model: (a) the spatial 
size of the neutrinos produced through a decaying particle with decay rate F is v/F for the 
parent particle at rest. This relation is commonly used in estimates of the neutrino wave 
packet sizes but this model provides an explicit example satisfying such relation, (b) The 
coherent creation of neutrino flavor states is possible because the kinematical contribution of 
neutrinos to their creation probabilities are negligible, allowing to define the flavor state as 
a superposition of mass eigenstates in the expected way. We should remark that the second 
conclusion is, in principle, applicable to more realistic decays because the main ingredient 
is universal: neutrino mass differences, compared to the charged lepton mass differences, 
contribute negligibly to their creation probabilities. 
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APPENDIX A: THE FUNCTION 



The function O appearing in Eq. (1261) can be written in the following simple form: 



% 3 



/ dv PJ VP> 1 iil (Al) 

2tt 2 ^W y E ] {E j + m J )M i + Ej- E„ v ; 



/OO ITT 2 

dE9{E - m 3 )^E 2 - m 2 (E - m 3 )f 2 {E) - j f- - (A2) 

7o^|^| 2 0i(^-M i ). (A3) 



7o ; 



In Eq. (1Aip we have assumed the cutoff function is isotropic /(p) = /(|p|), |p| = p. In 
addition, it is assumed in Eq. (1A2I) that / is a function of E only, i.e., /(p) — >■ f(Ej(p)). 
We used Eqs. (I3TJD. (|3T|) and 

C(P) = * a _ /( |^- =(-^)(p-^ 2 ) r , (A4) 
(2vr)2 x /2/; ;A //; ; • m ; 

^|^(p)| 2 = 2(/,(p)) 2 , (A5) 

where 

£(p) = 1 3 {fgL . (A6) 
If we extend the variables to small complex values, we obtain 



4>i(x±ie) =$ j (x)± l -f j (x), (A7) 



where 



fj(x) = 2n6(x - mj)(x - m^yj ' x 2 - m 2 f 2 (x) . (A8) 

all x while T'-{x) > for x > rrij am 
scale of /, assuming / is a smooth function on that region. 



We see Tj(x) > for all x while r^-(x) > for x > rrij and x much lower than the cutoff 



APPENDIX B: COMPLETENESS OF THE STATES \k(r), v 3 ;(p, s)) 



(±) 



We prove here the states \k(r), Vj(p, s)) (±) in Eq. ( 135|) complete the Hilbert space in the 



sector of one n or kiSj. We calculate 



ij ij 
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+ Yl J d3 ? rf Vl^,^'(q)>oA'j',u(q>q / )o(^,^(q / )l 

ij,i'j' 

+ J2 / rf 3 g^(q)|/ l ,z/ j (q)) 00 (7r| + h.c. (Bl) 

a 

The spin indices are suppressed. We get 

C = Z 2 J2 I ^I^(P)| 2 , (B2) 

Ai>j>,i»j»(q.,q[) = [^"^.^q^q) + a* Y/Y ,(q,q / ) + / d^a^fa, q)a^ iV ,(p, q')] (£3) 

i%(q) = Z*^(q) + W <fy4/$(pK,i'i'(P, q) • (B4) 
The first coefficient yields 

c = z 2 J rf 3 p|A,(p)l 2 



IVMi + ^p) +ie) 



= 7/ d£ £ 

Im /" 

A 



M£ + *e)| 2 



M 1+mi " \h (E + ie) 



1 



2?ri 
1 



Mi+mi 

dz 



1 



1 



/i (£ + «e) /io(f?-ie)J 



27ri J P ft<o(<z) 



= 1. 



(B5) 
(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 
(Bll) 



The contour P is a path along the real axis coming from 00 to m x + Mi below the real axis 
then going from mi + M 1 to 00 over the real axis. If the contour is closed with the aid of 
a very large circle, the integral over the closed contour is zero due to the absence of poles. 
The integral over the circle is equal to 1 due to ho(z) ~ z. 
We obtain ^^/^-//(q, q') = from 



/ ^«ij/i'(p>q)«^»i»(p,q') = d3p 



13 



X 



Mi + Ej(p) - M it - Ej/(q) + ie 
1 

Mi + Ej(p) - Min - E fl (q) - ie 



(B12) 



-1) 



dz 



2-ni J P h (z) z — M v — Ej>(q) + ie 
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Analogously one can prove 



z — Min — Ejn (q) — ie 

-9oVM)Vj(q!)Ui'j'Ut"j" 
{M v + E r (q) - M v , - E r (q>) - 
1 



ie 



h Q {M v + E r (q) - ie) h (M v , + E f ,(q) + 



a i'j',i"j 



»iq,q') - on"?',*? (q', q) 



(B13) 



(B14) 
(B15) 



S y (q) = 0. 



(B16) 



APPENDIX C: CALCULATION OF EQ. {77]) 

The gap between Eqs. (1761) and (1771) is filled by 

gUij f°° 4tt 9 /sin(pr)\ e -W+^-(p))t 



*r?( x ,t) = / dp< s (px) — — — ^ (Cl) 

J (2tt)I 7o ° V ; «9rV r J h(Ej + M t + ie) V 7 

4tt tft/y f°° dE rs . A , e"^' 

~ / n p,'X cos»,t- — ; — - — -, (C2) 

*r ( 27 r)f ^ ' J K 1 } 3 h(E) + fr(£) 1 J 



2v^ lr . f°° dE e~ iEt 

ir Jrrij+Mi V j E — + i^T it 

n fjl p -iM^t poo jp -iEt 

^U^yip^)-— dEcos((p J + -)r)-——. (C4) 



To get to Eq. (101 ft we made use of 

/^(p)^=r#^o*)^|(^). 



(C5) 



To get to Eq. (102)) we neglected the nonradiative term containing r 2 and used the change of 



variables p = |p| — > E = Ej(p)+Mi. We also used the shorthands pj = Pj(Ej) = y E 
vj = vj(Ej) = Pj( j^^ for Ej = E — Mi. To get to Eq. (103)1 we approximate 



m 



.i- 



h{E + ie)^E-M n + l -T v , (06) 

by Taylor expanding around E = M n . Notice that h(M n ) = and h'(M n ) = 1. The use of 
the approximation (106)) requires 

|r'(M 7r )|«i, |/ i "(M^)r 7r |«i. (07) 
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In Eq. 0C40 we approximated E = M w — Mj in the whole integrand, except in the expo- 
nentials and in the term h~ l . It is assumed that these functions vary slowly in the interval 
E pd M-n ± . For pj in the cosine we approximated 

Pj~Pj + ^, (C8) 

where E is obtained from the shift E —>■ E + M T compared to Eq. ( 1C3I) . It is also assumed 
that 

M n - T n > mj + M t , (C9) 

which allows us to replace the lower limit of the integral by — oo without changing the 
integral appreciably. At last, Eq. ( 1771) is recovered by using the integral 

-iEt 

dE - = -2m6{t)e- xt , A > , (CIO) 

j E + iX 

computed in the complex plane closing the contour on the upper or lower half plane. 
APPENDIX D: APPROXIMATE EIGENSTATE 



We can see the state in Eq. (1451) can not be an exact eigenstate of the total Hamiltonian H 
by calculating the eigenvalue equation, analogous to Eq. ( |25i) . with the complex eigenvalue 

E n = M n + i\. (Dl) 

We thus obtain 

M 7r + iA + o (M^ + zA) = 0. (D2) 
The real part is just Eq. (1251) . The imaginary part yields 

i r (M 7r ) i\ 
|A| -"2i + <W)--T' (D3) 

which can never be satisfied because of the minus sign. Therefore if both |A| and are 
small, the state (T43|) represents an approximate eigenstate of H. 

APPENDIX E: WIGNER-WEISSKOPF APPROXIMATION 



We can compare the results of Eq. (1551) with the transition probabilities with respect to 
the free (bare) states. Such transition amplitudes resemble the Wigner-Weisskopf (WW) 



approximation 



24[ | employed in atomic physics. 
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For the transition amplitude of the approximate pion state to the free pion we obtain 



Z 2 ~ 2 M<t)) = M0)|7T(t)> 



(El) 
(E2) 



where the factor Z 2 5 is necessary from renormalization. 

For the transition amplitude of the approximate pion state to the free l^Vj states we 
obtain 



MrlvfasMt)) = C3?(p) 



-i(Mi+Ej)t 



+ h(Mi + Ej + ie) / d V 



Cr. s '(p')\ 2 e~ i ^' +E ^ t 



M v + E v - Mi- Ej - ie 



-i{Mi+Ej)t 



1 _ e -i(M v -M^E 3 )t e -^ft 



E3) 
(E4) 



Equations (1E2I) and ( 1E4I) illustrate clearly the irreversible probability flow from the free state 
\tt)qZ^^ to the free states \k(r),Uj(p,s))o. 
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